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We describe nonlinear phonon-thermoelectric devices where charge current and electronic and 
phononic heat currents are coupled, driven by voltage and temperature biases, when phonon-assisted 
inelastic processes dominate the transport. Our thermoelectric transistors and rectifiers can be 
realized in a gate-tunable double quantum-dot system embedded in a nanowire which is realizable 
within current technology. The inelastic electron-phonon scattering processes are found to induce 
pronounced charge, heat, and cross rectification effects, as well as a thermal transistor effect that, 
remarkably, can appear in the present model even in the linear-response regime without relying on 
the onset of negative differential thermal conductance. 

PACS numbers: 73.23.-b,73.50.Fq,73.50.Lw,85.30.Pq 


I. INTRODUCTION 

Onsager’s formulation of irreversible thermodynamics 
and his famous reciprocal relations (1931)i, based on 
the time reversibility of microscopic dynamics, tie dif¬ 
ferent currents to driving forces (affinities) and provide 
strict bounds on the efficiency of energy transducers^. 
Among such transport phenomena the thermoelectric ef¬ 
fect which describes the coupling of electronic charge 
and heat currents nowadays attracts significant attention 
experimentally, computationally, and fundamentally, as 
new-unique thermoelectric materials and devices are pur¬ 
sued with a promise for an improved efficiency, see e.g. 
Refs. iS- The thermoelectric effect has been tradition¬ 
ally characterized by linear response quantities, Seebeck 
and Peltier coefficients. Nonlinear thermoelectric phe¬ 
nomena constitute a new area of research, anticipated to 
enhance thermoelectric responsei^rJ^. An elastic scatter¬ 
ing theory of nonlinear thermoelectric transport has been 
recently developed by considering self-consistent screen¬ 
ing potential o^^'^^~ — . Other studies had considered non¬ 
linear thermoelectric transport with explicit electron- 
electro n^^i^^ , electron-phononH and electron-photon^^ 
interactions. However, an investigation of nonlinear in¬ 
elastic “phonon-thermoelectric” systems, where charge 
current and electronic and phononic heat currents are 
coupled, and non-linear thermoelectric device operations 
(such as rectifiers and transistors) is still missing. 

By “phonon-thermoelectric” systems, we refer to a 
setup first considered in Ref. [13 (see also. Refs. IIB - IItI) . 
where electrical and thermal heat currents from source 
to drain were induced and manipulated by a third, in¬ 
dependent, phonon reservoir. This terminal is character¬ 
ized by its temperature, which is possibly distinct from 
the temperatures of the source and drain. The physi¬ 
cal mechanisms involving the phonon baths are inelastic- 
phonon assisted electron hopping from the source to the 


drain. In the setup shown in Fig. [TJa), the substrate 
corresponds to the third terminal, and it determines the 
temperature of phonons. In the linear response regime, 
a hot phonon bath can pump electrical current, to drive 
electrons against the electrochemical potential gradient. 
This effect was termed as the “three-terminal thermoelec¬ 
tric effect”, and it is similar to the photo-electric effect 
in photovoltaic systems where high-temperature photons 
pump electric currents. In fact, this analogy has been uti¬ 
lized for proposing useful thermoelectric devices^. Here, 
we extend this mechanism to the nonlinear response 
regime, to realize nontrivial thermoelectric functions and 
devices. 

Nonlinear electrical elements, diodes, amplifiers, and 
transistors, are key components in electronics, at the 
heart of modern technology. More recently, significant 
efforts have been devoted to the exploration of analogous 
phononic elements at the nanoscale^^. These endeav¬ 
ors, nonlinear nanoelectronic and nanophononic, were 
often separately pursued. However, it is to be noted 
that in many cases the interaction of electrons with the 
atomic degrees of freedom (vibrations, phonons) cannot 
be ignored. While electron-phonon dissipation effects of¬ 
ten impede device operation, here we exploit them to 
achieve compound nonlinear functionalities. We demon¬ 
strate these functionalities in a three-terminal geometry 
(the simplest situation). The effects described here can 
be generalized to other situations in which electronic cur¬ 
rents and phononic currents are spatially separated. 

In this work, we provide a unified description of cou¬ 
pled and nonlinear electrical and phononic transport 
in the three-terminal geometry when inelastic electron- 
phonon scatterings play a key role. The nonlinear, inelas¬ 
tic electron transport in our geometry (see Fig. [1} pro¬ 
vides charge and thermal rectification, cross-rectification 
effects (e.g., thermal rectification induced by voltage), 
as well as the realization of the thermal transistor effect 
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FIG. 1. (Color online) (a) Scheme of a DQD that can serve as 
a phonon-thermoelectric diode and transistor. The QDs are 
embedded in the nanowire and are controlled by gate volt¬ 
ages: I and r control the local potentials, t tunes the tunnel¬ 
ing between the QDs. The two electrodes, L and R, apply 
voltage and temperature biases across the QDs. The insula¬ 
tion layer suppresses the thermal contact between the metal 
electrodes and the substrate, which provides thermal energy 
to phonons, (b) Elastic and inelastic contributions to the 
charge current as a function of level detuning; microscopic 
processes are illustrated in panels (c) and (d). (c) Inelastic 
thermoelectric transport assisted by a phonon of frequency 
ujq. (d) Elastic thermoelectric transport. The shaded green 
area represents the broadening of the left and right QDs in¬ 
cluding the hybridization between them. The dashed arrows 
display the two main tunneling paths, of different energies. 
In (b) we used t = 15 /reV, = 1 MeV, ujq = 100 /reV, 
= 5 /reV, ksT = 20 fjeV (for Tl = Tr = Tph = T), 
Ee = 0, and hl = —fJ-R = 40 ^teV. 


and a phonon-thermoelectric cross effect (e.g., control¬ 
ling source-drain I^-V characteristics by modulating the 
temperature of a third terminal). Remarkably, we find 
that in the linear response regime inelasticity can lead to 
a thermal transistor effect. This is in contrast to common 
wisdo m^^i^^ that thermal transistor effect must comply 
with negative differential thermal conductance. The es¬ 
sential ingredient here is that inelastic electron-phonon 
scattering processes simultaneously involve three differ¬ 
ent reservoirs. This mechanism also works for inelastic 
electron transport assisted by coulomb drag effects^SiSi, 
magnons^^, photons^i^i, and plasmons^^, thus our re¬ 
sults are of interest beyond the specific model considered 
in this work. 


II. MODEL SYSTEM 

We consider a double quantum-dot (DQD) embedded 
in a nanowire in contact with two metals and a phonon 
substrate (“pli”), see Fig. [TJ The quantum dots (QDs) 
are defined by voltage gates (labeled by I and r in Fig. 
HJa) with tunable electronic energy levels Ei and Ep. t 
is a hopping element between the QDs and are the 
hybridization energies of the dots to the source and drain 
electrodes, labeled by L and i?, respectively. Charge cur¬ 


rent, electronic heat current, and phononic heat current 
are induced by applying a voltage bias between the ter¬ 
minals L and R and a temperature difference between 
the three terminals. This setup (and related models) 
have been explored in e.g., Refs. 123 - 1^ . How¬ 

ever, previous studies had focused on the linear response 
regime, while here we uncover coupled nonlinear phe¬ 
nomena. The system is described by the Hamiltonian, 

EI — HrjqR) -\- Re—ph T Eigad Htunei -b Hpfi 

with 

Hdqd = ^ Eiclci -b (tcjcr + H.C.), 

i—^,r 

^ ^ UqC^Cfi^CLq H.C., 

Hlead = 

j—L,R k 

Hfunel — E VL,kc\cL,k + y]] VR^kc\.CR^k + H.C., 

k k 

Here c] {i = £, r) creates an electron in the i-th QD with 
an energy Ei. The £ (r) QD is located next to the left 
(right) lead. The tunneling elements from the £ QD to the 
right lead and that from the r QD to the left lead are as¬ 
sumed negligible. Electron-phonon interactions (matrix 
element Uq) allow inelastic electron transport through 
the system; aj creates a phonon with a wavevector q 
and frequency ujq. We adopt the conventions that h=\ 
and ks = I throughout the paper. To maintain a finite 
temperature difference between the electrodes and the 
substrate, we suggest to thermally isolate them with a 
layer of thermal and electrical insulation [see Fig. [ija)]. 
The phonons involved in the inelastic transport through 
the DQD system may be confined phonons in the quan¬ 
tum wire, or bulk phonons in the substrate. We assume 
the thermal contact between the quantum wire and the 
substrate to be good so that these two types of phonons 
acquire the same temperature Tph. The phonon tem¬ 
perature can be controlled by the substrate, while the 
temperature of the electrodes (Tr and Tr) can be indi¬ 
vidually controlled by heating them with, e.g., AC elec¬ 
tric fieldsii^. Although Coulomb interaction is certainly 
important in electron transport effects through QDs, to 
simplify the problem and to uncover the essential physics, 
we assume that the intra-dot Coulomb interaction energy 
is much larger than other relevant energy scales in our 
problem, as well as that the inter-dot Coulomb interac¬ 
tion is small and negligible. Under these assumptions, we 
can ignore the Coulomb interaction altogether. This as¬ 
sumption limits our discussions to the regime where there 
is no more than one electron in each QD. Nevertheless, 
the essential element for the phenomena described in this 
work is energy exchange between electrons and phonons, 
rather than electron-electron energy exchange process. 
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FIG. 2. (Color online) Charge, heat, and cross rectification 
effects, (a) Charge current Is for Ee = —Er = —2fcsr and (b) 
optimization of charge rectification Re with QDs energies Ee 
and Er at Sfi — 20 /reV. (c) Electronic heat current Iq (”aW” 
denoting atto-Watts) for Ee = —bksT and Er = —ksT, and 
(d) optimization of the associated rectification strength Rt, 
by tuning Ee and Er using 5T = O.ST. (e) Electronic heat 
current as a function of applied voltage for Ee = —ksT and 
Er = 2fcsT, and (f) optimization of the relevant rectification 
Rte at the bias 5/r = 20 ^teV. (g) Charge current against the 
(metals) temperature difference for Ee = —^ksT and Er = 
—ksT, and (h) optimization of the rectification strength Ret 
by controlling the DQDs at ST = O.ST. We used t = 15 fieV, 
$0 = 1 fJ-eV, LOO ~ 100 /reV, je = Jr = 5 /xeV, ksT = 20 /xeV, 
and Tph = To for all figures. 


III. CURRENTS IN AND BEYOND 
LINEAR-RESPONSE 


Non-equilibrium steady state quantities of interest are 
the electric current Ig , the electronic heat current travers¬ 
ing from the left lead to the right lead Iq = ^{—Ql + 


Qr), and the phonon heat current I^ = —Qph, with Qi 
{i = L,R,ph) denoting the heat current flowing into the 
*th reservoir. Both elastic and inelastic processes con¬ 
tribute to transport in the system [see panels (b), (c) 
and (d) in FigUj. 

The inelastic (“mef) contribution to the currents is 
calculated from the Fermi Golden-Rule treatment of Ref. 
setting for convenience p = + tr) = 0, 

^e\inel — \inel — ■ ( 3 ) 


Here e < 0 is the charge of an electron, E = 
n = Er — Ee. In the above equation, /jv = ^ 


with Te^r = 'ye-phfd^-fr)Np and F^-^^ = je-phfr{i- 
fe)N^. Here = Nr + ^ ± 4sgn(H) with the Bose- 
Einstein distribution for phonons Nr = [exp(^P^) — 1]“^. 


For the commonly used zinc blend semiconductors such 
as GaAs and InP, the electron-phonon interaction en- 


,.33.46-48 


ergy Uq determines the transition rat o'^'^i"'-’ — je-ph = 


^o(!f)"exp 


. Here stands for the electron- 


phonon scattering strength, n provides the power-law de¬ 
pendence on phonon energy H with a characteristic en¬ 
ergy loq. Appendix A provides details on the deriva¬ 
tion of this expression from the microscopic theory of 
electron-phonon interactions in GaAs QDs. It is found 
that in this system the piezoelectric mechanism domi¬ 
nates the electron-phonon interaction. In this situation, 
we adopt the parameters <^0 = 1 and n = 1—4^. 
The Fermi Golden Rule method is valid when kRT ^ 
IItIt ^ iF^-^r — These conditions can be main¬ 

tained by tuning the coupling between the QDs and the 
leads. In this regime, the steady state distributions on 
the two QDs can be approximated by those of the nearby 
leads2^, i.e., fe ~ fL{Ee) = ) -I- 1]“^ and 

fr ^ fR[Er) = [exp(-^i^^) -b I]-\ where /l and Jr 
are the Fermi distribution functions for the left and right 
electronic reservoirs, respectively. If the Fermi Golden 
Rule assumptions (specifically, the second inequality) are 
not satisfied, the inelastic current can be calculated more 
generally via a rate equation technique [see Appendix B]. 
The Fermi Golden Rule method is exploited below for 
the analytic study of transport coefficients. Numerical 
results in this work are calculated directly from the rate 
equation method (unless specified otherwise). 

The contribution of elastic (“el”) processes to the cur¬ 
rents is given by Landauer’s formula 


f dr 

Ie\ei = e J ^T{e)[fL{e) - /^(e)], (4a) 

lQ\ei = J ^ene)[h{e)-fR{e)], I^q^i = 0. (4b) 


Note that elastic processes do not contribute to the 
heat current from the phonon terminal I^. The energy- 
dependent transmission function is obtained from the 
Garoli formula, T(e) = tr[G''(e)fj;,G“(e)f^j, where 
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G’-(e) = G“t(£:) = 


Tl = 


It 0 
0 0 


e - + -t 

— t e — Er + 

0 0 ■ 


-1 


, and Ffl; = 


0 7^. 


where we assume 


constant (energy-independent) tunneling rates and 
7 r. This results in T(e) = t^ 7 ^ 7 r/|dGp with do = 
{e — Ef, + i^i/2){e — Er + i^r/‘^)—t^■ In Fig. [T](b) we plot 
the elastic contribution to the charge current (green) as a 
function of the DQD level detuning for an otherwise sym¬ 
metric junction 7 j = jr- As expected, the current is an 
even function with respect to detuning. In contrast, the 
inelastic contribution (blue) is asymmetric with respect 
to detuning. This distinction serves as a signature of in¬ 
elastic transport in DQD systems in experiments^^. Con¬ 
tributions from elastic transport processes are included 
in our simulations. 

The entropy production rate for the whole system is 
given by 


_ Ql Qr Qph 




(5) 


where the affinities for /e, /q, and [identihed in 
Eq. (ini) by Ii,l 2 ,h separately], are Ai = (^ + 

^l_'i A _ _J_ _ J_ A — 1 I 1 _ 1 .^P- 

2TRh ^2 — Tr Tl^ ^3 — 2Tl ^ 2Tr Tp^ ’ 

spectively. By setting the reference temperature at Tq = 

Tl+Tr ’ obtain Ai = ^, A 2 = XrTr ’ A 3 = TqTp 

where 5yL = fiL-fJ-R, ST = Tl-Tr, and AT = Tph-Tc^. 
Note that the reference temperature is defined according 
to ^ = i(^ -I- ^) [differing from the average temper¬ 
ature T = ^{Tl + Tr) in nonlinear-response regime]^. 

Expanding Eqs. © and (|4]) to capture nonlinear ef¬ 
fects, we find, to the lowest nontrivial order 


li — MijAj -h Lijf^AjAj^ + 0{A^), (6) 

j jk 


where My = Mf^ + M™®* with My denoting the linear- 
response coefficients. The Onsager reciprocal relations, 
My = Mji, hold for both elastic and inelastic transport 
processes. The second-order terms Tyfe only come up 
from the inelastic transport processes. Elastic coefficients 
are given by 

Mli = ^j dene)r‘>{e)[l - r‘>{e)], (7a) 

Mil = M 2 I = ^J dene)r{e)[l - r{e)]e, (7b) 

M 2 I = ^J dene)r‘>ie)[l - rie)]e\ (7c) 

M|'=M;'=0, j = 1,2,3. (7d) 


The superscript eq denotes the equilibrium distribution. 
The linear response inelastic coefficients which satisfy the 
Onsager reciprocal relations are 

M™^' = s.s.rZr, ( 8 ) 


where T^^p = “ fr‘^)iNp is the tran¬ 

sition rate from QD I to QD r at equilibrium, and 
Si = e, S 2 = T, S 3 = D. 

The second-order transport coefficients are given by 

where (gjk = gkj) 

911 = - f!% 

gi2 = eT(/;« - /D + ieO(/r + - 1), 

5l3 = 0, 323 = 0, 

322 = (E^ + io2)(/-^ - /;«) + TD(/;« + /;«-1), 

333 = 2 D|D |(1 + 2 A^«). ( 10 ) 

Note that Onsager reciprocal symmetry does not guar¬ 
antee Tijk — Tji}^^. 

For the described functionalities to be prominent, the 
contribution from the inelastic processes needs to be en¬ 
hanced while the elastic processes must be suppressed. 
This is realized by strong electron-phonon interaction, 
high temperature^^, and the mismatch of QDs energy 
levels. The latter creates a barrier for elastic tunnel¬ 
ing and favors inelastic transport. We noticed from Fig. 
mb) that there is an optimal alignment of QDs energy 
levels for inelastic processes to be dominant. It is also 
noted that the inelastic contribution is nonzero even at 
D = 0 , which is a special property of the piezoelectric 
electron-phonon interaction; The electron-phonon cou¬ 
pling strength is proportional to D, but this linear depen¬ 
dence cancels out the divergence in the phonon number 
in the limit D —>■ 0 , overall leaving a finite term which 
contributes to the current. A careful examination of the 
D —0 limit for realistic systems is beyond the scope 
of this paper. Our main conclusions do not rely on the 
behavior of the system in this special regime. 

The transport coefficients depend on the average en¬ 
ergy E (representing deviation from particle-hole sym¬ 
metry) and the detuning D (measuring deviation from 
mirror symmetry). The symmetry of transport coeffi¬ 
cients and currents under (i) particle-hole transforma¬ 
tion: Tph '■ E —>■ —E (hence E —>■ —E and D —D) and 
(ii) mirror symmetry (assuming 7 ^ = 7r) 7 m '■ Ec Er 
(i.e., D —D) are summarized in Table I. These sym¬ 
metries rule out certain nonlinear effects in the presence 
of mirror symmetry and/or particle-hole symmetry. For 
example, when E = 0 but D 7 ^ 0 only terms of the same 
symmetry under the particle-hole and mirror transforma¬ 
tions remain nonzero. 


IV. RECTIFICATION (DIODE) EFFECT 

Coupled thermal and electrical transport allows un¬ 
conventional rectihcation, for example, charge rectifi¬ 
cation induced by a temperature difference, besides 
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TABLE 1. Symmetries of Qij, Lijk, Mij, and the heat currents. 


Terms Particle-hole Tph Mirror 7 m 


Pll, 922 , 933 

Odd 

Odd 

512 

Even 

Odd 

Till, L 122 , Li33, 7/212 

Odd 

Odd 

7/112, 7/211, 7/222, 7/233 

Even 

Odd 

7/311, 7/322, 7/333 

Even 

Even 

7/312, 7q^, Mi2 

Odd 

Even 

CO 

Odd 

Odd 

M 23 

Even 

Odd 

Mil, M 22 , M 33 

Even 

Even 


TABLE II. Functionality of second-order coefficients 


Terms (Lijk) 

Diode or Transistor effect 

7/111 

charge rectification 

7/222, L333 

electronic and phononic heat rectification 

7/233, 7/322 

off-diagonal heat rectification 

7/122, Li33 

charge rectification by ST or AT 

7/211, 7/311 

heat rectification by voltage 

7/212, 7/112 

other nonlinear thermoelectric 

7/321 

phonon-thermoelectric transistor 


showing the standard charge and heat rectification ef¬ 
fects. The magnitude of the rectification effect is de¬ 
fined by Re = charge rectification, 

Rt = II+ -ST )I (electronic) heat rectifica¬ 
tion, Ret = I+[/ ^-Jt )I charge rectification in¬ 

duced by the temperature difference dT, and Rte = 
|/e'^(v)|+|j^e (-v)l foi' heat rectification induced by voltage 
bias. We shall refer to Ret and Rte as “thermoelectric 
rectifications”. The dependence of these quantities on 
the DQDs energies is displayed in Fig. [21 demonstrating 
significant rectification effects. 


From these figures, one finds that thermal diode Rt 
and thermoelectric diode Rte cannot be realized when 
the system is invariant under Tph ® Tm (i-e., vanishing 
E). In Table II we classify second-order transport co¬ 
efficients as rectifiers and transistors. We note that the 
rectification and transistor behavior as displayed in Fig|2] 
[obtained from a direct calculation using Eqs. ©-(j!])], 
acquire the same symmetry as the corresponding second- 
order transport coefficients listed in Table llll Specifically, 
Re and Ret (relating to Lm and L122, respectively) are 
odd under both Tph and Tm, whereas Rt and Rte (relat¬ 
ing to L222 and T211, respectively) are even under Tph, 
but odd under Tm- In addition, it is demonstrated in 
Fig|2] that one can tune En and E^ for reaching optimal 
electrical, thermal, and cross rectifications. 



FIG. 3. (Color online) (a) Linear thermal transistor. We plot 
the electronic heat currents Iq and Iq (the former describes 
the current leaving the L terminal) and the phonon heat cur¬ 
rent (”aW” denoting atto-Watts), against the tempera¬ 
ture of the phonon bath. The two electronic heat currents 
vary significantly whereas the phonon current changes very 
little with Tph- Ei = A.^ksT, Er = SksT, ST = 0.02r, 
kpT = 20 /xeV, and dp = 0. (b) Phonon-transistor. The 
charge current T, plotted as a function of dp, is largely con¬ 
trolled by the temperature of the phonon bath. Ei = —AksT, 
Er = 2kBT, and ST = 0. Other parameters are the same as 
in FigE] 


V. TRANSISTORS 

One of our central results is that a thermal transistor 
effect can develop in the linear-response regime. Specifi¬ 
cally, the heat current amplification factor is given by 



M23 — 5M33 

Ee 


M33 

n 


where Iq = Iq — /q^/ 2 is the heat current flowing out of 
the L electrode (since p^ = p^). Eq. m follows directly 
from (i) Iq = EilH + ‘'‘'elastic contributions'', I^ = HIn 
[according to Eq. ©], and (ii) the elastic contribution 
to Iq does not depend on the phonon temperature Tpt- 
Hence 

dTphlq = EtdppjN, = Q-dT^jN, ( 12 ) 

and one obtains Eq. m- Remarkably, heat current am¬ 
plification, a > 1, is achievable once \Eg\ > |n| in both 
the linear and nonlinear regimes whenever the inelastic 
contribution is nonzero. This conclusion counters present 
designs in which a thermal transistor effect can only be 
realized in the nonlinear response regime based on a neg¬ 
ative differential thermal conductance^^. In fact, one can 
prove that a thermal transistor effect can materialize in 
the linear response regime only when inelastic processes 
involve more than two reservoirs due to restrictions im¬ 
posed by the second law of thermodynamics [see Ap¬ 
pendix cp-. The inelastic electron hopping mechanism 
assumed in our three-terminal setup is one of the sim¬ 
plest examples which can realize such a nontrivial prop¬ 
erty: electrons enter from the source and leave at the 
drain with the assistance of phonons from a third termi¬ 
nal (the substrate). Fig©a) demonstrates the thermal 
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transistor behavior: Iq changes more significantly than 

Specifically, for the adopted temperature difference 
ST = 0.02ro ^ To and for small AT = Tph — To the 
system is indeed described by linear response terms. Al¬ 
though the coefficient a increases for smaller fl, the vari¬ 
ation of the heat current Iq as controlled by the temper¬ 
ature of the phonon terminal can also decrease once the 
inelastic transport mechanism weakens [see Fig. 1(b)]. 

In addition to the thermal transistor effect exemplified 
in Figl31[a), we also show in Figl31[b) that by changing 
the temperature of the phonon bath we can substantially 
modify the I^-V characteristics — This effect is analo¬ 
gous to the field effect transistor, but here it is controlled 
by the phonon temperature instead of the gate voltage 
at the third terminal. This functionality is also differ¬ 
ent from the unconventional transistor effect observed in 
superconductor-normal-metal-superconductor junctions 
by Saira efc al. where gate-voltage at the third ter¬ 
minal controlled thermal transport between source and 
drain using the Coulomb blockade in the normal-metal 
region^. 


here, i.e., the essential importance of inelastic transport 
effects for a linear-response transistor and other nonlin¬ 
ear device functionalities, is also of importance for related 
research fields. 
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VI. CONCLUSIONS AND DISCUSSIONS 

We proposed a realistic and relatively simple setup for 
the realization of thermal, electrical, and cross rectifiers 
and transistors by exploiting phonon-assisted hopping 
transport in DQD systems in a three-terminal geometry. 
While we used a DQD system embedded in a phonon sub¬ 
strate as our platform, results are applicable to other two- 
site (two-level) systems maintained out of equilibrium by 
multiple biases and subjected to inelastic transport pro¬ 
cesses. Numerical simulations using realistic parameters 
demonstrate strong charge, heat and thermoelectric rec¬ 
tifications. Furthermore, for the first time, it is observed 
that a thermal transistor can be realized without negative 
differential thermal conductance. In addition, a transis¬ 
tor effect where the le-V characteristics is tuned by the 
phonon temperature was revealed. These functionalities 
should enable smart manipulations of heat and charge 
transport in nano-systems, core elements in future infor¬ 
mation processing technologies. Although calculations 
were performed for a specific set of parameters, the un¬ 
covered functionalities should be observed whenever in¬ 
elastic transport processes dominate elastic effects. As 
discussed in Refs. [H and [H, this can be achieved in 
the high temperature and strong electron-phonon inter¬ 
action regime when elastic transport is suppressed e.g. 
by a barrier, or an energy gap. Our study here essen¬ 
tially aimed at the hybridization of two distinct branches 
of technology: electronics and phononics^^. Amalgama¬ 
tion of technologies may provide new platforms and op¬ 
portunities for high-performance, high-energy-efficiency 
nanotechnologies. Extensions to consider the interplay of 
spin, charge, and heat transport in the presence of inelas¬ 
tic effects are of interest for exploring spin-caloritronics^ 
devices beyond linear response. The physics revealed 


Appendix A: Derivation of the phonon induced 
transition rate between DQDs. 


The electron-phonon coupling matrix for longitudinal 
phonons in GaAs is given by 




2 pMVsiqV 


- iq^), 


(Al) 


where V is the volume of the system, pM = 5.3 x 
10^ kg/m^ is the mass density of GaAs, S = 7.0 eV 
is the deformation potential, and Vsi = 5.29 x 10^ m/s 
is the velocity of longitudinal phonons. The piezoelectric 
potential is 


T^gT = 


87rei4 {iqxqyqz) 


q^ 


(A2) 


where ei 4 = 1.41 x 10® V/m is the piezoelectric constant 
and At = 12.9 is the relative dielectric constant. Those 
material parameters are adopted from the standard semi¬ 
conductor handbook. Ref. l53l. From the above equations, 
one finds that 


ia,d^ = 


32It7r®e^ej4 {hqxqyqzY 

K^PMVsiqV g® 




2pMVslV 


(A3) 


When considering the coupling of electrons to transverse 
phonons, only the piezoelectric mechanism contributes. 
In this case the two transverse branches of phonons give 
the total contribution 


ICufP = 


32;t7r"e^ej4 

K^PMVstVq 


qiqt + qiqi + qtqi 


i^qxqyqz 


(A4) 
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where Vgt = 2.48 x 10^ m/s is the velocity of transverse 
acoustic phonons. 

We now consider electronic wavefunctions in the quan¬ 
tum dots, assuming a Gaussian form, 


4-1 (f) = 


4-2 (r) = 




(27r)3/4/3/2 


exp 


exp 


+ y2 [z + d/2f 


472 


4Pj 


+ 2/2 {z-d/ 2 y 


4/2 


4/2 

^^qd 


where Iqd is the characteristic length of the wavefunction. 
The two wavefunctions are centered at (0,0, ±(i/2) (i.e., 
the distance between the two QDs is d), respectively. 



FIG. 4. (Color online) (a) Electron-phonon coupling strength 
vs. quantum dot radius Iqd for longitudinal phonon 
with deformation potential coupling mechanism (“LDP”), 
transverse phonons with piezoelectric coupling mechanism 
(“TPE”), and longitudinal phonons with piezoelectric cou¬ 
pling mechanism (“LPE”). Calculated from Eqs. (E3-(Il9l) 
using material parameters of GaAs. (b) Characteristic 
phonon frequency ujx as functions of quantum dot radius Iqd 
(“L” and “T” stand for longitudinal and transverse phonons, 
respectively). 


With these elements at hand, we write down the tran¬ 
sition rate between the dots using the Fermi Golden Rule, 


qX 

X (5(±/7a7A,5 ~ Ej- + El). 


The integral over space gives 
|(4'i|e*'J'''|4'2)P = exp 


■ ^ 4 


(A5) 


(A6) 


for transverse (longitudinal) phonons with the piezoelec¬ 
tric coupling mechanism, ojldp = ^lpe = Vsi/lqd, 
UJTPE = Vst/lqd, riLDP = 3, riLPE = UTPE = 1, and 


^LDP 


ihPE 


irPE 


2'KhpMvllld ( 4 / 2 ^) ’ 

967re2ef4 / \ 

i^hK’^PMV^lqd y 4/2^ j ’ 

1287re2e24 / \ 

ihhp'^PMvltlqd J ' 


(A7) 

(AS) 

(A9) 


We calculate for GaAs QDs using d/lqd = 4. Results 
are plotted in Fig. SI For Iqd — 12 nm, ^tpe = 1 peV and 
lot PE — 120 /teV, approving the parameters, E^nd ioq 
(representing ^tpe and lotpe in the main text, since we 
consider only the piezoelectric electron-phonon coupling 
for transverse phonons there), chosen in the main text. 


Appendix B: Comparison between the Fermi-Golden 
Rule approximation and the Rate equation method 


The Fermi Golden Rule approximation is valid when 
ksT ^ 7^,7r ^ iFf^r — Tr^e^. Here we will examine 
situations in which the condition ^ ~ 

is relaxed. We assume that the broadenings of the QD 
levels 7^,7r are still much smaller than the thermal en¬ 
ergy ksT. The transport currents can then be calcu¬ 
lated through Eq. ©• However, the particle current In 
itself attains a more complex form in terms of the elec¬ 
trochemical potentials and temperatures. In this regime, 
the steady state distributions on the two QDs, fg and fr, 
can be calculated by solving the following rate equations 
in steady state. 


0 = -^ = -idfe - fL{Eg)] - -fe-ph 

X [/,(! - fr)N- - Ml - fg)N+l (Bl) 

0 = ^ = -7r[/r - fpiEr)] + le-ph 

X [Ml - fr)N- - /,(! - fg)N+]. (B2) 

Once obtaining fg and fr, the particle current can be 
calculated via 

In = le-ph[fM - fr)N- - MI - MN+]. (B3) 


The summation over q can be converted into integration 
which gives 


^e—ph 


E 

X=LDP,TPE.LPE 



Here LDP stands for longitudinal phonons and the de¬ 
formation potential mechanism, and TPE {LPE) stands 


When 7 ^, 7 t. ^ \In\, the Fermi Golden Rule approxima¬ 
tion is validated, since fg ~ fpiEg) and fr ~ fpiEr). 

We plot le-V characteristics for different phonon tem¬ 
peratures using the rate equation method and the Fermi 
Golden Rule approximation in Fig. [5] It is noted that 
differences are negligible unless we operate in the far- 
from equilibrium regime, Tp^ ^ Tt,Tn, or ^ kgT. 
Our results in Fig. [5] indicate that the nonlinear perfor¬ 
mance is marginally reduced; qualitative features remain 























FIG. 5. (Color online) I^-V curves for different phonon tem¬ 
peratures calculated from the rate equation method (denoted 
as “RATE” in the figure) and the Fermi Golden Rule (denoted 
as “FGR”) approximation. The parameters are the same as 
in Fig.EIb). 


the same. We speculate that the phenomena described in 
this work persist beyond the strict regimes where the rate 
equation method or the Fermi Golden Rule approxima¬ 
tion are justifiable. Simulations in the main text of this 
paper were calculated using the rate equation method. 


Appendix C: Restrictions on the thermal transistor 
effect in the linear-response regime from the 
second-law of thermodynamics 


Here we discuss restrictions, directly arising from the 
second-law of thermodynamics, on the thermal transistor 
effect which we materialize in the linear-response regime 
[see Fig. EKa)]. Our discussion below is based on gen¬ 
eral irreversible thermodynamic arguments, and is holds 
for any thermodynamic systems in the linear response 
regime. Since our analysis here is confined to the thermal 
transistor effect, we shall consider pure thermal conduc¬ 
tion (i.e., the electrochemical potential difference is set 
to zero). 

For a system with three reservoirs, one can define three 
heat currents (one for each reservoir). However, due to 
energy conservation, Iq+ Iq += 0, we are left with 
two independent heat currents^i^. If we choose the heat 
current flowing out of the source and that emerging from 
the phonon terminal as the two independent currents, 
the thermal transport equation for three terminal in the 
linear response regime can be generally written as. 



(Cl) 


Here Iq = Iq — I^ /2 is the heat current flowing out 
of the source. T is the equilibrium temperature. 
and Kph are the diagonal thermal conductance, while Ko 
stands for the off-diagonal thermal conductance. The 
above equation is derived using Onsager’s theory of 
linear-response by finding the thermodynamic forces con¬ 
jugated to the two heat currents. It can also be obtained 
by linearizing the transport equation in the main text. 


Eq. ([6]), and then reorganizing the transport coefficients 
according to the definition of Iq. The second-law of ther¬ 
modynamics imposes the following restrictions on trans¬ 
port coefficients, 

Kl > 0, Kph > 0, KhKph > Kl (C2) 


We now note that in a three-terminal system, one can 
generally identify two distinct classes of transport mecha¬ 
nisms: (i) those involving only two terminals (reservoirs) 
in each microscopic process (briefed as “2T-class”); (ii) 
those in which energy is exchanged among three reser¬ 
voirs in each transport process (briefed as “3T-class”). 

For the 2T-class, the thermal transport equations in 
the linear response regime are generally given by, 

Il^r = Klr{Tl - Tr)/T, (C3a) 

iL^ph = KLphiTL - Tph)/T, (C3b) 

iR^ph = KRph{TR — Tph)/T, (C3c) 

where Krr, K^ph, and KRph are the thermal conduc¬ 
tances between pairs of reservoirs. The above equations 
represent Fourier’s law of thermal conduction. The sec¬ 
ond law of thermodynamics requires that the thermal 
conductances are positive, 

Klr>Q, KLph>0, KRph>0. (C4) 

We can organize equations (IC3|) into the general form 
m, written in terms of a new set of transport coeffi¬ 
cients, 

K'l = Klr + Klph, 

K'ph = KLph + Krph, K = -KLph- (C5) 

Based on the restrictions (IC4I) . we find that [on top of 
Eq. (IC^ j. 

K<0, (C6) 


As discussed in the main text, the heat current am¬ 
plification factor in the linear response regime is given 

by 






K 


ph 


(C7) 


From Eq. (IC6I) it is now clear that there is no thermal 
transistor effect for 2T-class models in the linear-response 
regime since they identically result in a < I. (Recall, the 
2T-class transport mechanisms correspond to a three ter¬ 
minal situation, with transport processes occurring inde¬ 
pendently between every pair of terminals. ) In contrast, 
for 3T-class transport mechanisms, the restriction (IC2I1 
does not disallow heat current amplification, a > 1. In 
fact, as discovered in this work, one can manifest the 
thermal transistor effect in the linear response regime for 
the phonon-assisted hopping process, which is an exam¬ 
ple of the 3T-class transport mechanisms. 
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